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ABSTRACT 

We determine the effect of rotation on the luminosity of supermassive stars. We apply the Roche 
model to calculate analytically the emitted radiation from a uniformly rotating, radiation-dominated su- 
permassive configuration. We find that the luminosity at maximum rotation, when mass at the equator 
orbits at the Kepler period, is reduced by ~ 36% below the usual Eddington luminosity from the corre- 
sponding nonrotating star. A supermassive star is believed to evolve in a quasistationary manner along 
such a maximally rotating "mass-shedding" sequence before reaching the point of dynamical instability; 
hence this reduced luminosity determines the evolutionary timescale. Our result therefore implies that 
the lifetime of a supermassive star prior to dynamical collapse is ~ 56% longer than the value typically 
estimated by employing the usual Eddington luminosity. 
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1. INTRODUCTION 

Recent observations provide strong evidence that super- 
massive black holes (SMBHs) exist and are the sources 
that power active galactic nuclei and quasars (see, e.g., 
Rees 1998, for a review and references). However, the sce- 
nario by which SMBHs form is still very uncertain (see, 
e.g., Rees 1984, for an overview). One promising route is 
the collapse of a supermassive star (SMS). Once they form 
out of primordial gas, sufficiently massive stars will evolve 
in a quasistationary manner via radiative cooling, slowly 
contracting until reaching the point of onset of relativis- 
tic radial instability. At this point, such stars undergo 
catastrophic collapse on a dynamical timescale, possibly 
leading to the formation of a SMBH (Bisnovatyi-Kogan, 
Zel'dovich & Novikov 1966; Zel'dovich & Novikov 1971; 
Shapiro & Teukolsky 1983). 

Because most objects formed in nature have some an- 
gular momentum, rotation is likely to play a significant 
role in the quasistationary evolution, as well as the final 
collapse of a SMS. The slow contraction of even a slowly 
rotating SMS will likely spin it up to the mass-shedding 
limit, because such stars are so centrally condensed. At 
the mass-shedding limit, matter on the equator moves in 
a Keplerian orbit about the star, supported against grav- 
ity by centrifugal force and not by an outward pressure 
gradient. The SMS evolves in a quasistationary manner 
along the mass-shedding curve, simultaneously emitting 
radiation, matter and angular momentum until reaching 
the onset of radial instability. 

In this paper, we derive the luminosity of a uniformly 
rotating SMS as a function of its spin rate, up to the mass- 
shedding limit. The magnitude of the luminosity is cru- 
cial because it determines the evolutionary timescale of the 
star as it evolves. Elsewhere we use the result to follow the 
slow contraction of a cooling, rotating SMS to the onset 
of dynamical instability (Baumgarte & Shapiro 1999, here- 
after Paper II). Here, however, we focus on the emitted flux 
and total integrated luminosity from a stationary SMS as 
a function of its rate of rotation. [For an overview of pre- 
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vious work on SMSs and references, see, e.g., Zel'dovich & 
Novikov 1971; Shapiro & Teukolsky 1983; and Paper II.] 

Even though our calculation, which is analytic up to 
quadrature, is rather simple and straight-forward, we have 
not been able to find a similar argument in the literature. 
Previous analytical arguments have dealt with more gen- 
eral rotation laws, but they adopt the slow rotation ap- 
proximation and emphasize gas-pressure atmospheres (see, 
e.g., Kippenhahn 1977). While detailed numerical calcula- 
tions of luminosities of rotating stars have been carried out 
for select main sequence and massive stars (see, eg, Tassoul 
1978, Table 12.1, and references therein; Langer & Heger 
1997), we cannot find a calculation for a SMS. Hence in- 
dependent of its relevance to the evolution of SMSs prior 
to catastrophic collapse, our result may be of interest to 
stellar modeling of rapidly rotating stars in the limit of 
very high mass, where our calculation is applicable. 

Our paper is organized as follows: in Section 2 we enu- 
merate and justify our assumptions. In Section 3 we briefly 
review the Roche approximation, which we use to describe 
the outer layers of a rotating SMS. In Section 4 we derive 
the flux and luminosity from the star. In Section 5 we 
briefly summarize our results and compare them with pre- 
vious calculations of rotating main sequence stars. 

2. BASIC ASSUMPTIONS 

Our analysis relies on several explicit assumptions, all 
of which we expect to hold to high accuracy in SMSs. In 
particular, we assume that the star is 

1. dominated by thermal radiation pressure; 

2. fully convective; 

3. uniformly rotating; 

4. characterized by a Rosseland mean opacity that is 
independent of density; 

5. governed by Newtonian gravitation; 
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6. described by the Roche model in the outer envelope. 

For large masses, the ratio between radiation pressure, 
P r , and gas pressure, P g , satisfies 
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(see, e.g., Shapiro & Teukolsky 1983, eqs. (17.2.8) 
and (17.3.5)); here the coefficient has been evaluated for 
a composition of pure ionized hydrogen. For stars with 
M ;> 10 4 M Q , we can therefore neglect the pressure con- 
tributions of the plasma in determining the equilibrium 
profile, even though the plasma may be important for de- 
termining the stability of the star (Zel'dovich & Novikov 
1971; and Shapiro & Teukolsky 1983). A simple proof that 
SMSs are convective in this limit is given in Loeb & Rasio 
(1994). This result implies that the photon entropy per 
baryon, 
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is constant throughout the star, and so therefore is j3 » 
8(s r /fe) _1 . Here a is the radiation density constant, ns is 
the baryon density, and k is Boltzmann's constant. As a 
consequence, the equation of state of a SMS is that of an 
n = 3 polytrope: 



P = Kp A ' 3 , K = 
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where P is the pressure, p the mass density, m the atomic 
mass unit and p the mean molecular weight (cf. Clayton 
1983, eq. 2-289; note that Clayton adopts a different defini- 
tion of @, which is related to ours by /^ciayton = 

The third assumption, that the star is uniformly rotat- 
ing, is probably the most uncertain of our assumptions. 
Nevertheless, it has been argued that convection and mag- 
netic fields provide an effective turbulent viscosity which 
dampens differential rotation and brings the star into uni- 
form rotation (Bisnovatyi-Kogan, Zel'dovich & Novikov 
1967; Wagoner 1969). 

In the high temperature, low density, strongly ionized 
plasma of a SMS, Thomson scattering off free electrons is 
the dominant source of opacity. This opacity is indepen- 
dent of density and justifies our fourth assumption. 

We assume that gravitational fields are sufficiently weak 
so that we can apply Newtonian gravity. SMSs of interest 
here have R/M ;> 400 (see Paper II), so that this assump- 
tion certainly holds. Relativistic corrections are important 
for the stability of SMSs, but can be neglected in the anal- 
ysis of the equilibrium state. 

Finally, the Roche approximation provides a very ac- 
curate description of the envelope of a rotating stellar 
model with a soft equation of state, as in the case of an 
n = 3 polytrope (for numerical demonstrations, see, eg, 
Papaloizou & Whelan 1973 and Paper II). Since our anal- 
ysis is based on this approximation, we will briefly review 
it together with some of its predictions in the following 
section. 

In applications to SMSs our analysis neglects electron- 
positron pairs and Klein-Nishina corrections to the elec- 
tron scattering opacity, which is valid for M ;> 10 5 M Q 
(see, e.g., Fuller, Woosley & Weaver, 1986). 



3. REVIEW OF THE ROCHE MODEL 

Stars with soft equations of state are extremely centrally 
condensed: they have an extended, low density envelope, 
while the bulk at the mass is concentrated in the core. 
For an n — 3 polytrope, for example, the ratio between 
central density to average density is p c /p = 54.2. The 
gravitational force in the envelope is therefore dominated 
by the massive core, and it is thus legitimate to neglect the 
self-gravity of the envelope. In the equation of hydrostatic 
equilibrium, 



VP 

= -V($ + $ c ), 

p 
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this neglect amounts to approximating the Newtonian po- 
tential $ by 

M 

d> = 5 

r 

(here we adopt gravitational units by setting G = 1). 
In ([|) we introduce the centrifugal potential <& c , which, 
for constant angular velocity f2 about the z-axis, can be 
written 



-ftV sin 2 (9. 
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Integrating eq. (Q) yields the Bernoulli integral 

h + $ + $ c = H, 
where H is a constant of integration and 



dP , P 
— = n + 1 - 
P P 



(6) 
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is the enthalpy per unit mass. Evaluating eq. (R) at the 
pole yields 



Rn 



(9) 



since h — on the surface of the star and $ c = along the 
axis of rotation. In the following we assume that the polar 
radius R p of a rotating star is always the same as in the 
nonrotating case. This assumption has been shown nu- 
merically to be very accurate (e.g., Papaloizou & Whelan 
1973). 

A rotating star reaches mass shedding when the equator 
orbits with the Kepler frequency. Using eqs. (||) and (0), 
it is easy to show that at this point the ratio between 
equatorial and polar radius is 



(10) 



L P/ shedd * 

The corresponding maximum orbital velocity is 
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(Zel'dovich & Novikov 1971; Shapiro & Teukolsky 1983). 
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Fig. 1. — The flux F as a function of polar angle for different values of a = Sl/f! s h e dd- Note that at the mass-shedding limit a = 1 the flux 
vanishes on the equator. 



4. LUMINOSITY OF ROTATING STARS 

According to our assumptions, the pressure in super- 
massive stars is dominated by radiation pressure 
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In the diffusion approximation, the radiation flux is every- 
where given by 

F = -^-VU. (13) 
Here U the energy density of the radiation, 



U = aT 4 



3P, 



(14) 



and k is the opacity (which we assume to be dominated 
by electron scattering, n = k cs ). Inserting ( |l_3| ) and (14) 
into the equation of hydrostatic equilibrium Jm yields 



kF = V($ + $ c ). 



(15) 



In polar coordinates in an orthonormal basis, the magni- 
tude of the flux is 



F = [F~ 
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Evaluating the gradients of $ and $ c in the envelope yields 

21 V2 

• (17) 



F = 
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Introducing the dimensionless spin and radius parameters 

n 



hcdd 



and 



and denoting 
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the usual Eddington flux from a spherical star, we can 
rewrite equation (Il7j) as 
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Note that, from eq. (Q), the surface of the star is defined 

by 

$ + $ c - H = 0, (22) 

or equivalently 2 
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(23) 



Given a and 8, the value of z on the surface can be found 



by solving this cubic equation. Eqs. (15) and (E3) imme- 



diately imply that the flux is normal to the surface of the 
star. 

Evaluating eq. (pTl) at the surface, we plot the emer- 
gent flux F as a function of 6 for different values of a 
in Figure 1. Note that at the mass-shedding limit, when 
a = 1, the flux vanishes at the equator (where z — 3/2 and 
sin# = 1). This, of course, is an immediate consequence of 
hydrostatic equilibrium: at mass shedding, the centrifugal 
force exactly balances the gravitational force at the equa- 
tor, so that the pressure gradient vanishes (eq. (11)). For 
radiation-dominated stars, eq. ( p"5| ) then implies that the 
flux has to vanish. 

The total luminosity can be found by integrating 



L 



F dA 



A 



FdA 



(24) 



over the surface A of the star. The surface element dA 
can be written 



dA = 2irr sin Ods = 2irr sin (dr 



l d9 u 



>l/2 



/ \ dr 
= 2irr 2 sin 6d9 1 + — — 
\ r z dO 

or 

dA = 2nr 2 dn (l + (1 - , 

Here we have introduced 

/i = cos 9 



1/2 



1/2 



and 



z 



dz 
d\i 



(25) 



(26) 



(27) 



(28) 



FEdd 

2 Expanding eqs. ( |l7j ) and (E3h to lowest order in a 2 shows that they are in perfect agreement with eqs. (39) and (30) in Kippenhahn (1977) 
for uniform rotation. Kippenhahn's treatment allows for non-uniform rotation, but is restricted to slow rotation. 
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Fig. 2. — The luminosity L as a function of the orbital velocity Q. At the mass shedding limit, the luminosity is reduced by 36 %. 



Differentiating (p3[), the latter can be expressed as 

8 z 3 a 2 n 



z 
z 
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Putting the pieces together, we find that the luminosity is 
given by 



2 J 2nr z dn ( 1 



'\ 2 
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or 
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where I/Edd is the usual Eddington luminosity 
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It proves most convenient to evaluate eq. (31) numerically. 
In Figure 2, we plot the resulting luminosity L as a func- 
tion of the spin parameter a. Obviously, for nonrotating 
stars with a — we recover L = L^dd- For maximally ro- 
tating stars, however, the luminosity is reduced by about 
36%: 

-^shedd = 0.639-LEdd- (33) 

Accordingly, adopting the Eddington luminosity for a su- 
permassive star that evolves along the mass shedding limit 
would underestimate its lifetime by about 36%. 

5. DISCUSSION 

We find that the luminosity of a SMS rotating at break- 
up velocity is reduced by about 36% from the luminosity 
of a nonrotating SMS of the same mass. 



It is difficult to compare this result with previous nu- 
merical calculations of massive, rotating stars, which are 
summarized in Table 12.1 in Tassoul (1978, compare the 
discussion in Kippenhahn 1977). No calculations seem 
to have been performed for stellar masses greater than 
62.7M Q . The luminosities of these 62.7M Q models at 
break-up velocity are indeed reduced below the nonrotat- 
ing, spherical luminosities, but only by 7%, much less than 
what we find. However, the physical conditions in 62.7M Q 
stars are very different from those in SMSs and do not sat- 
isfy our assumptions (see Section 2). For example, at these 
moderate masses, the stars are not dominated by radiation 
pressure; according to eq. (|), j8 » 1 for these stars and 
P varies with both the location in the star and the orbital 
velocity 3 . Also, the total opacity, at least close to the sur- 
face of the moderate mass stars considered previously, is no 
longer dominated by electron scattering and contains non- 
negligible contributions from bound-bound and bound-free 
absorption. These contributions introduce a dependence 
on density, so that our assumption 4 no longer holds. More 
specifically, spinning up the star may decrease the density 
in the envelope, therefore decrease the opacity, and hence 
increase the luminosity. This would partly compensate for 
the decrease of the luminosity due to rotation, and will 
decrease the effect that wc find for true SMSs. Finally, we 
note that moderate mass stars are not fully convective. 

Our result is important because the luminosity deter- 
mines the timescale of the evolution and hence the life- 
time of SMSs, which are believed to evolve along a mass- 
shedding sequence. This lifetime has been used in several 
calculations of SMS evolution. We adopt our new result in 
Paper II, where we analyze the secular evolution of SMSs 
up to the onset on radial instability. 
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